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Top-down algorithm 
based on MODIS Fire 
Radiative Power 
(AQUA/TERRA) 


FRP Emission factors 
tuned by means of 
inverse calculation 
based on MODIS AOD 
data. 


Daily mean emissions, 
NRT 


Prescribed diurnal 
cycle 

In GEOS-5 BB 
emissions are 
deposited in the PBL. 


Aerosol Data Assimilation 


Q State representation QO Observation operators 
> Multiple 3D concentrations > Intrinsic aerosol optical properties needed for 
» Mass remotely sensed data: 
» Number (modal schemes) » Mass extinction coefficient, single scattering 


» Bin sizes (Sectional schemes) aoeOO Diase anal 


+ Niumibeiror traearenteneio bundieds » These are often poorly known but assumed 


to be known due to identifiability issues: 


T=BeM 


O Emissions: 


» Dynamic: dust, marine, biogenic aerosols 
> Given diversity of aerosol representation in 
» Remotely sensed: biomass burning 
models, DA software should not hardwire 


» Inventories: anthropogenic optical properties. 


Aerosol Observing System 


QO Aerosol Optical Depth (AOD) is the 


most commonly available observable | arc) (AR ee ae 

> Vertically integrated mass weighted by aa | (Observing system + > 
extinction coefficient, summed over = ae 7 ce05-5 nerasois : 
multiple species: low observability eee z a 


> Available multi-spectral AOD 
measurements are not really measured 
Q Radiance assimilation: 
> Vector scattering calculations needed for O Lidar measurements provide vertical info 
UV-VIS measurements are not cheap 


> Surface BRDF characterization is a 
challenge > Attenuated backscatter again requires optical 


> Spatially coverage is poor (pencil thin) 


QO Surface PM 2.5 assumptions which are not directly measured 
> Single level 
> Often plagued by representativeness 


» New HSRL concept is promising 
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2D AOD ANALYSIS 


Observable 550 nm AOD is 2D 


Constrains column averaged optics 


Cannot constrain speciation or vertical 


distribution 


Analysis in observation space: 


T* = Hq" =H (q’ + 4q") 
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GOING TO 3D CONCENTRATIONS 
Based on error covariances: 
5q* = BH? (HBH*) 574 
Using ensemble perturbations, 
dq? = XYT (YYT) 579 
NRT GEOS-5 uses Local Displacement 
Ensembles (LDE), in aD 


gq” MMR Update 
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ORIGINAL MODIS C6 AOD BIAS CORRECTED AOD 
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MODIS Neural Net AOD Retrievals trained on AERONET 


MODIS Collection 6: Aqua 
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APPLICATIONS 
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Unique amongst its peers, the MERRA-2 reanalysis now includes an 
aerosol reanalysis for the modern satellite era (1980 — onward). 
Aerosols are coupled to the meteorological reanalysis (both radiatively 
and through emissions/loss processes). 


MERRA-2 Global, Monthly Mean Aerosol Optical Depth (AOD) 
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First aerosol assimilation to include major historic volcanic events like El Chichon 
(1982) and Pinatubo (June, 1991). 

Movie shows the co-evolution of gaseous SO2 emissions from Pinatubo (left) and 
formation of the the sulfate aerosol plume (right) as SO2 is converted into particles. 
SO2 (g) is from emissions inventories and unconstrained by assimilation. Sulfate 
aerosol AOD (right), however, is impacted by the assimilation of total aerosol AOD. 


Field Campaign 
Support 


Global 5-day chemical - 

forecasts , \w 
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field missions including TC4, 
ARCTAS , GloPac, ATTREX, Time (Hours UTC) 
DISCOVER-AQ, HS3, 
iS) oN GVA, ts¥ 40) OS NOF 
ORACLES, etc. 


22 


Comparison of observed (top) and simulated (bottom) 
aerosol backscatter for a flight during the 2013 
SEAC4RS campaign. 
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Aerosol Impact on Meteorological Assimilation 


QO While hybrid aerosol-meteorology data 


assimilation may not be computationally 
attractive in the near future, separately 
analyzed aerosol fields may be used as 
parameters in the observation operator 
of IR channels (AIRS, IASI, etc.) 


QO Figure on the right shows monthly mean 


brightness temperature difference 
(aerosol minus no-aerosol) for several 
GSI IR sensors with fractional dust AOD 
> 0.65. 


See Jong Kim poster for more details. 
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Mean A Tb, Du Frac > 0.65, 20160801-31 (12z) 


Mean A Tb, Du Frac > 0.65, 20160801-31 (12z) 


—1.0 -0.5 
A Tb (°K: Aer-Ctl) 


, SPACE FLIGHT CENTER 


GOING FORWARD 
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Aerosol EnKF fas 


QO As part of GMAQO’s hybrid system, 
aerosol ensemble members are 
produced as a matter of routine 


O The same Whitaker-Hamill EnKF used ° 
for the hybrid Meteorological 
assimilation has been adapted for 


PSAS AOD inc (current method) EnKF AOD inc Ensmean 


PSAS - 20161216 - 9z EnkF (2D) - 20161216 - 9z 
a a A ors = 


aerosols 
O Target observation systems , | | : : | | ,: 
> Multi-spectral AOD: 470, 550 and 870 nm POP eg at nae ae de km yeuae es Mae carga nee 


> Lidar attenuated backscatter 


> Sensors: MODIS, VIIRS, GEO, CATS/CALIOP, 
TropOMI 1 dared 
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Concluding Remarks 


QO Aerosols are an integral part of the GEOS-5 N.R.T. and re-analysis systems 
>» Aerosols are run at the same resolution as the main meteorological model 


O MERRA-2 provides the first integrated aerosol-meteorology reanalysis for the satellite 
era 


QO Current GEOS-5 developments incorporate cloud and aerosol microphysics 

>» Aerosol-cloud interactions, better representation of aerosol mixing state 

> There is great need for in-situ and remotely sensed data to evaluate/validate new parameterizations 
0 QFED biomass burning emissions to include diurnal cycle from GEO satellites 
QO Aerosol analysis migrating to an EnKF based system 

> New observables: multi-spectral AOD, attenuated backscatter 


> New active/passive sensors: VIIRS, GEO, CATS/CALIOP, TropOMI 
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